Modified Scaling Relation for the Random-Field Ising Model 
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We investigate the low-temperature critical behavior of the 
three dimensional random-field Ising ferromagnet. By a scal- 
ing analysis we find that in the limit of temperature T — > 
the usual scaling relations have to be modified as far as the 
exponent a of the specific heat is concerned. At zero temper- 
ature, the Rushbrooke equation is modified to a + 2/3 + 7 = 1, 
an equation which we expect to be valid also for other systems 
with similar critical behavior. We test the scaling theory nu- 
merically for the three dimensional random field Ising system 
with Gaussian probability distribution of the random fields 
by a combination of calculations of exact ground states with 
an integer optimization algorithm and Monte Carlo methods. 
By a finite size scaling analysis we calculate the critical expo- 
nents v M 1.0, f3 « 0.05, 7 « 2.9 7 w 1.5 and a » -0.55. 



Above two dimensions, the ferromagnetic random-field 
Ising model has an ordered phase for low temperatures 
and small random-fields as was proven by Imbrie Q and 
later also by Bricmont and Kupiainen For larger 
fields the system develops a domain state J3| which has 
been shown to have a complex and fractal structure [Q. 
It is now widely believed that the phase transition from 
the ordered to the disordered phase is of second order. 
In three dimensions, the values of some of the critical 
exponents are now well established, like /3 w 0.05 and 
v « 1. Although real space renormalization yields de- 
viating results concerning v (see e. g. ||) the values for 
(3/v are in the same range. However, a complete set of 
values of the critical exponents fulfilling the predicted set 
of scaling relations could not be established, neither 
by experimental measurements (for a review, see ||) per- 
formed usually on diluted antiferromagnets in a magnetic 
fields which are thought to be in the same universality 
class nor by any numerical methods (see e. g. |pH). 
Especially, the value of a - and even its sign - is highly 
controversial. 

The Hamiltonian of the RFIM in units of the nearest 
neighbor coupling constant J is 



H 



ai<T i - H B * 



(i) 



The first sum is over the nearest neighbors and the 
spin variables <Tj are ±1. The random- fields Bi are 
taken from a Gaussian probability distribution P(Bi) ~ 
exp(-(B l /A) 2 /2). 

We assume that there is a zero temperature fixed point 
at a finite value A r of the random-field width. Introduc- 



ing the scaling variable / = A — A — g(T) (see also fg]) 
where the condition / — describes the critical line we 
expect the same critical behavior no matter if we vary the 
temperature or the random-field strength. Hence, for the 
singular part of the internal energy E it should be 
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and for the singular part of the free energy F 
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In the limit of low temperatures F equals E and the 
question arises which derivative - first or second - with 
respect to the random-field strength yields the exponent 
of the specific heat? 

The scaling ansatz for the singular part of the free 
energy for zero homogenous magnetic field is: 



F(T,f) = \f\ 1/x2 T d 



T 



(4) 



Hence, for the most singular part of the specific heat 
it is c ~ l/l 1 ^ 2-2 and consequently a = 2 — l/x 2 , as 
usual. On the other hand, for fixed critical random-field 
A = Ao and in the limit of temperature T — ► the 
prefactor T in Eq. || becomes critical and hence it is a — 
1 — 1 /%2- This is also consistent with the scaling behavior 
of II = -^{F + TS) with S = -ff which follows from 
the scaling ansatz for the free energy. The most relevant 
singular terms are: 
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In the limit T — > the first term vanishes and only the 
second term is observed in the specific heat leading to a. = 
1 — \/x2 as above. Consequently, it follows by standard 
scaling theory that the scaling relations in the limit of 
temperature T — > have to be modified with respect to 
a. For zero temperature, the equation corresponding to 
the Rushbrooke equation has the form 



2/3 + 7 = 1 



(6) 



Additionally, it is remarkable that in Eq. || the sec- 
ond more singular term is small for either small T or 
small d\f\/dT. Therefore, one can expect to observe the 
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anomalous zero temperature critical behavior as long as 
the critical line is flat. 

A similar, although more complicated consideration 
holds for the specific heat. Building the derivative c = 
-T§y§, the most relevant singular terms are: 



for the magnetization and 
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The most relevant term has the square of the slope -g^ 
as a prefactor while the next relevant term has the cur- 



vature 



a-l/l 

dT 2 



as a prefactor. If the critical line starts hor- 
izontal at T — but with a finite curvature the most 
relevant term will be suppressed and the unusual less 
critical behavior will be observed, yielding a = 1 — l/x2 
for low temperatures. Note that in both cases, for c as 
> critical behavior is an inflection 



well as for §£ the T 
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point since and ^{ change the sign at the critical 
point. However, for finite temperatures close to the crit- 
ical point a crossover to the "normal" critical behavior 
can be expected. 

In order to test these arguments numerically we con- 
sider the three dimensional RFIM and calculate exact 
ground states (EGS) using an optimization algorithm 
well known in graph-theory. The Ising-system is mapped 
on an equivalent transport network, and the maximum 
flow is calculated using the Ford-Fulkerson algorithm 
fL2| |l4|| . We used a simple cubic lattice with periodi- 
cal boundary conditions and linear lattice sizes varying 
from L — 6 to L = 20. From the spin configurations 
of the ground state, we can calculate the magnetization 
M = [m] av , where m = j^J2i a ii the internal energy, 
E = [/i] av where h = j-sH and the disconnected suscep- 
tibility Xdi s = L 3 [w 2 ] ay , where the square brackets de- 
note an average taken over 30-1900 random-field configu- 
rations, depending on the system size. The advantage of 
the numerical technique above is that it supplies equilib- 
rium information. But on the other hand it is restricted 
to zero temperature. Therefore, we combine the EGS cal- 
culation with Monte Carlo (MC) methods. Starting at 
zero temperature with an EGS-spin configuration for a 
certain set of random-fields we heat the system slowly us- 
ing the standard heat bath algorithm. Hence, we get MC 
data at low temperatures, T < T c /2, which are close to 
equilibrium. We checked that by heating the system with 
decreasing heating rates until no further change of the 
data was visible. Using this MC method we can addition- 
ally calculate the susceptibility x = ^p" [(m 2 ) — (to) 21 

and the specific heat c = |4 [(h 2 ) — (h) 2 
angles denote a thermal average. 

From the scaling relations above follow the finite size 
scaling relations 



I av 



av 
where the 



M = L- f3 ' y M ((A - AJL 1 ^) 



(8) 



Xdis = L^X ((A - A C )L^) 



(9) 



for the disconnected susceptibility. 

Figure [l] shows the scaling plot for the magnetization 
data from EGS calculations as described above yielding 
A = A C (T = 0) = 2.37 ± 0.05, v = 1.0 ± 0.1, and f3 = 
0.05±0.05. These are values which are not surprising and 
in agreement with most of the previous work, especially 
the previous EGS calculations of Ogielski | TJ] . The error- 
bars are estimated since there is no straight-forward way 
to extract error-bars from a finite-size scaling plot. 
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FIG. 1. Scaling plot of the magnetisation from EGS. 

Using the same values for A c and v as above from 
the scaling plot for the disconnected susceptibility (not 
shown) we get 7 = 2.9 ± 0.3 which is also in agreement 
with most of the previous work. 

The first derivative of the ground state energy shows 
the critical behavior of the specific heat. The behavior of 
E can be understood by a series expansion of the energy 
in the vicinity of the critical point: 



E(A) =E a + E 1 {A- A c ) + E S (A - AJ 1 " 



(10) 



The i?i-term is important since as we argued above a can 
be expected to be negative for low temperatures following 
Eq. |(| Hence, the finite size scaling form is 
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for the derivative of E in the critical region. Differentiat- 
ing our energy data numerically, we obtained the scaling 
plot shown in Figure ^. Once more we used the same 
values for A c and v as above and chose E\ such that 
|j2| = E\ at the inflection point. Note, that we neglect 
here a possible size dependence of the analytic parts of 
the energy, i. e. a possible L-dependence of E\ which 
is obviously very small as Figure || demonstrates. This 
analysis leads to a — —0.55 ± 0.2. 
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FIG. 2. Scaling plot of the derivative of energy from EGS. 

The approach to determine a for finite temperatures 
is a direct MC simulation of the specific heat. Figure || 
shows the corresponding data. The ground state spin 
configurations were used as initial spin configurations 
for a MC simulation - of course of systems with iden- 
tical random-field configuration. Then the systems were 
slowly heated (10000 MCS per temperature with tem- 
perature steps of 0.2). Data are shown for T = 1.4 
which is roughly 30% of the critical temperature at zero 
field. We do not find any divergence of the specific heat, 
i. e. no size dependence of the maximum of c. Hence, 
as above we analyzed the data subtracting the value of 
the energy at the inflection point cq. We took the values 
A C (T = 1.4) = 2.35 and v = 1 from MC simulation data 
of M for the same temperature (data not shown here). 
Our analysis yields once more a = — 0.55±0.2. The value 
A C (T — 1.4) is very close to the zero-temperature value 
Ao, confirming that the critical line is nearly horizontal 
in the low temperature region. Hence, as discussed above 
the true critical behavior fulfilling the Rushbrooke equa- 
tion is hard to observe and within our numerical accuracy 
we can only find the zero temperature exponent. 
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FIG. 3. Scaling plot of the specific heat from MC simula- 
tion. T = 1.4. 



]l5| ]. We analyzed MC data also for other temperatures 
but we did not find any significant temperature depen- 
dence of the critical exponents in the temperature range 
< T < 2 which is nearly one half of the phase diagram. 

To summarize, the values we determined for the crit- 
ical exponents of the 3D RFIM are in good agreement 
with many of the previous works, experimental as well as 
theoretical. The exponents 7 and 7 fulfill the Schwartz- 
Soffer equation 7 = 27 [ p6[ . The modified hyperscaling- 
relation p7| which can be written in the form (without 
a) 7 = Dp — 2/3 is also fulfilled by our exponents. 

It is the aim of this work to calculate as many ex- 
ponents as possible independently in order to test Eq.^ 
which is the most important aspect of our work. It is de- 
rived from the standard scaling ansatz for the zero tem- 
perature fixed point. Therefore we expect this equation 
to be valid also for other systems for which the discussed 
scaling ansatz is true. Candidates may be the random- 
field- Heisenberg model in appropriate dimensions. For 
the RFIM in higher dimensions new results suggest that 
there is a break of universality (i. e. the critical behavior 
depends on the kind of the distribution of the random- 
fields) as was shown in ref. [ jl8[ for four dimensions and 
earlier in refs. and p3] for the mean field solution 
of the RFIM. Additionally, it was shown [^lj that there 
is replica symmetry breaking for the mean field solution 
of a random field model with m-component-spins in the 
limit of large m. However, the replica symmetric solu- 
tion of the RFIM with a Gaussian distribution of random 
fields has a zero temperature fixed point and indeed, in 
the limit T — > the mean-field exponents fulfill eq. ^ 
since it is (3 = 0.5, 7 = 1 and a — —1 which can be in- 
ferred from a remark in ref. p9[ | stating that the entropy 
vanishes as T — > linearly with T" . 

We argued that the crossover to the normal critical 
behavior might be hard to observe as long as the criti- 
cal line is horizontal. We directly determined the con- 
troversially discussed exponent a for zero temperature 
yielding a — —0.55 confirming the validity of EqJ^. The 
same value is also observed for finite but low tempera- 
tures, the crossover to the true critical behavior cannot 
be observed within our numerical accuracy. Surprisingly, 
this value is also in agreement with recent Monte Carlo 
simulations p"T[ as well as with recent experimental re- 
sults p^] . Both did not find a divergence of the specific 
heat but a < 0, although these measurements and sim- 
ulation, respectively, were performed for higher temper- 
atures where true critical behavior should be easier to 
observe. 
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A standard finite size analysis of the MC data for x 
yields 7 = 1.5 ± 0.2 (not shown). This value is in agree- 
ment with the previous work of Ogiclski. Nevertheless 
there are deviations from previous results from series 
expansion where 7 was found to be significantly higher 
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